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We consider HQET including the first order correction in 1/w. A strategy for the computation of 
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is discussed. Only two quantities <f>i/2 have to be considered in order to match QCD and HQET, 
since the spin-dependent interaction is easily eliminated due to the spin symmetry of the static 
theory. Quite simple formulae relate the renormalization group invariant b-quark mass (Mj,) to the 
B-meson mass. All entries in these formulae are non-perturbatively defined and can be computed 
in the continuum limit of the lattice regularized theory. For the numerically most critical part, we 
illustrate the cancellation of power divergences by a numerical example. 
Numerical results for the 1/m correction to Mb, are presented in a companion talk. 

XXIIIrd International Symposium on Lattice Field Theory 

25-30 July 2005 

Trinity College, Dublin, Ireland 



*This work has been supported by the Deutsche Forschungsgemeinschaft (DFG) in SFB Transregio 9 "Computa- 
tional Particle Physics". 
' Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



A strategy for the computation of nth including 1 /m terms 



Rainer Sommer 



1. Introduction 

Although HQET is the most natural effective theory for heavy-light systems, its lattice reg- 
ularized version has practically only been used at lowest order. Indeed, a strategy to overcome 
the problem of power divergent mixings [|l|], was only found rather recently iQ]. Its potential was 
demstrated by a computation of the b-quark mass to lowest non-trivial order in 1 / m, the static 
approximation. Here we fill the formalism of [^, sketched in the abstract, with practicable defini- 
tions in terms of Schrodinger functional correlation functions and give a concrete formula for the 
1 / m-correction to the quark mass. 

Neglecting 1 /m^ corrections - as throughout this report - we write the HQET Lagrangian 

^HQET = -Sf'stat {x) - ft>spin ^spin {x) - Ofen ^km {x) ( 1 ■ 1 ) 

^spin = ^aBlZ/h , ^kin = ^DVh , (1-2) 

such that the classical values for the coefficients are (0\aa = Wspin = l/(2m). Since expectation 
values 

{ff) = + (,^„(^)kin + Wspin(^)spin , (1-3) 

(^)kin = £(^^kin W)stat , (^)spm = £(^^spin(x))stat (1-4) 

-Y X 

are defined through insertions of the higher dimensional terms i^kin, i^spin in the static theory, they 
are renormalizable by power counting. However, in order to have a well defined continuum limit 
the bare, dimensionful, couplings cOkim ^spm have to be determined non-perturbatively [[T], ^. In the 
framework of lattice QCD, this is possible by matching a number of observables, O,, / = 
between QCD and HQET, thus retaining the predicitivity of QCD. It is essential to note that this 
matching can be carried out in a finite volume of linear extent Li ~ 0.4fm, where heavy quarks can 
be simulated with a relativistic action [Q, ||, Q]. 

Since the lowest order theory is spin-symmetric, it is trivial to form spin-averages which are 
independent of Wspin- One thus expects that « = 2 is sufficient for a computation of the quark mass 
(in addition to (0\aa there is an overall (state-independent) shift of energy levels, which we denote 
by nibare)- For unexplained notation we refer to 

2. Basic observables 

We consider the spin-symmetric combination 

/r(0,r) =Zj{/i(75)}>/^{/i(7i)}3/\ (2.1) 
formed from the boundary to boundary correlation functions 

= L (Ci'(u)rcUv)?b(y)r?i(z)) , (2.2) 

u,v,y,z ^ ' 

of the QCD Schrodinger functional of size T xLr' and a periodicity phase d [§| for the quark fields. 
Replacing the b-quark field by the effective field y/h, using eq.( |L3|jL4| ), and accounting for the 
multiplicative renormalization of the boundary quark fields , one finds the \/m expansion 

/r = Z2 Z^e-'--^ {/r + o^nff""] , (2.3) 
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where the aformentioned energy shift mbare enters. Deviating from the choice in [Q], we now define' 

<I>i(L,M) = ln(/r(e,r)//r(0',r)) -\n{ff-\d,T)/fr\e' ,T)) (2.5) 

for r = L/2, 

^2{L,M) = j-ln{fr{e,T-a)/f^^e,T + a)) , (2.6) 
with the expansion 

^1 (L,M) = w^Rf- , ^2{L,M) = L (mbare + + CO^nTf'') (2.7) 

„k:n_ /i'""(e,r) ff-{e\T) 

rr = iiln{fr{d,T-a)/fr\d,T + a)) , (2.9) 

^\ff^\d,T-a) ff\d,T + a))- ^ ' 



T^kin 
A 1 



3. Step scaling functions 

We choose Li « 0.4fm, where a computation of <I>;(Li,Mb) is possible in lattice QCD (while 
at significantly larger values, Li/a would have to be too large in order to control effects). From 



eq. ( |2.7| ) one then gets ft)kin5"^bare for lattice spacings a = ^ x 0.4fm. On the other hand, contact 
to physical observables, e.g. the B-meson mass is made in large volume, where finite size effects 
are exponentially small. For reasonable values a/L\ = 1/12 and Loo ~ 1.5fm at the same lattice 
spacing, one needs L„/a ~ 50. This situation is avoided by first computing step scaling functions 
which connect <I>,(Li,M) to <I>;(L2,M),L2 = 2L\ and then connecting to large volume. 

With the Schrodinger functional coupling, u = g^{L), everywhere, the continuum step scaling 
functions a are defined by 

7?'f"(2L) 



(2L,M) = ar («)<I>i (L,M) , af" («) = Hm 

and 



(3.1) 



„=g2(L) 



^2(2L,M) - 2<I>2(L,M) = Omiu) + [o^„2L(rf"(2L) -rf"(L))] (3.2) 

rf"(2L)-rf"(L) 



am(M) + a|^'"(M)<I>i(L,M), o^"{u) = Urn 2L 



a/L^O /?ifi"(L) 
Here the static step scaling function 



u=fiL) 



a^{u) = \\m2L [rf (2L) - TT\L)\ ^^^..^^ , (3.3) 

is not identical to Omiu) defined earlier [Q], since Fj'^' differs from F**^^' defined there. Note that 
the step scaling functions are independent of M, but <I>,(L,M) have a mass dependence from fixing 
<I>,(Li ,M) in the full theory. 

' In the static computation of |^ tlie logaritlimic derivative T of tiie correlation function /a of the axial current with 
a boundary operator was used as a quantity to match effective theory and QCD. Including 1 /m terms its expansion reads 

f _ yHQETy y „-"ibi„Ao / rftat , HQET rftat , „ j-kin , „ rSpin"! r-^ a\ 

,/a - ''fi/f'^ |/a +'-A /sA+^in/A + fspin/A ]■ (2-4) 

with the term /J™ due to the 1/m correction to the static axial current. While CUspin represents no problem, an extra 
observable is needed to fix /^'a • Here, we avoid this complication by working exclusively with /j"*. 
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4. Large volume 

The connection of <I>, to the spin-averaged B -meson mass, m^, is 
LmB - <I>2(L,M) = [L{E'^' -Tf\L))] + [L(0^,{E''- -T^\L))] 



(4.1) 



[L{E''^'-Tf\L))\+p{u)^^{L,M), p{u)= lim L 



7Jdn T-kin 



rf"(L) 



Here we have used the abbreviations 



limrf (L), 



^kin 



lim rr(L) 



/?f"(L) 



u=f(L) 



(4.2) 



where E^^^^ is the (unrenormalized) energy in large volume in the spin-averaged B -channel in static 
approximation and ft)kin£''^" is its \/m correction. The hat on f"^'" is to remind us that this quantity 
turns into an energy only upon multiplication with the dimensionful (0\aa- Its numerical evaluation 
has already been investigated in [0]. We use [...] braces to indicate combinations which have a 



continuum limit by themselves. For example, the two terms in eq. (4.1) can be computed with 
different regularizations if this is useful. 



5. Final equation 



The above equations are now easily combined to yield the 1/m correction, to the (spin- 
averaged) B -meson mass via (L2 = 2Li), 



Stat , (Ifl) , (l^') 



(5.1) 

L2mf\M) = [L2 (E^'^' - rr (L2))] + aUui) + 2<I>2(Li,M) (5.2) 
L2m'i''\M) = a^''{m)^i{LuM), Ui=f{Li) (5.3) 
L2m'i'\M) = [L2(£"^"-rf"(L2))o;kin] =p(«2)al^"(«i)<J>i(Z.i,M). 



Again, terms in braces have a continuum limit. While tn^"^ is purely derived from finite volume, 
the term m^^^ involves a large volume computation. 

Starting from M^^^\ the solution of the leading order equation, 



and the slope 



d 

dM 



S = ^mf' 



Id, 



(1) 



(5.4) 
(5.5) 



we finally obtain the first order correction to the RGI b-quark mass 



1 



-m 



(5.6) 



The final uncertainty for Mb due to the 1/m expansion is of order 0(Aq(^p/M^), which translates 
into a numerical estimate of MeV scale. It is thus clear that other sources of error will dominate in a 
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practical calculation. Note that the precise value for me matters. One should use the spin-averaged 
mass ^^''P'^'-™^"'^' = i^go + 3^g, = [^5279+ |5325]MeV = 5314MeV if one can extrapolate E 
to the chiral limit of the light quark or 

^experimental ^ + |mB5 - fmBo = [5370 + | ( 5325 - 5279)] McV = 5405 MeV (5.7) 

if one works directly with a strange quark (as light quark). The latter formula neglects the depen- 
dence of the spin splitting on the light quark mass. 

6. Remarks 

The following facts are worth noting. 

• The 1 /m expansion in heavy light systems is an expansion in terms of Aqcd / where all 
external scales have to be of order Aqcd- This applies in particular to our scale L^^ . Indeed, 
numerically it is rather close to Aqcd and explicit investigations [||, ^] have shown that the 
1 /ni-expansion is well behaved even when L^Ms a factor two larger. 

• In our static computation [Q, we made the more natural choice T instead of Fi. Although 
it is advantageous to use Fi when one includes the \/m terms, the strategy can easily be 
formulated with F, at the expense of introducing a third quantity <I>, to fix c^*^^^. Since this 
will certainly be required for the computation of the 1 /m-correction to Fb, we will follow 
also that approach. 

• Note that at each order k in the expansion, the result is ambiguous by terms of order \/m^^^ . 
Thus both M^'' and M^'^' have an order \/m ambiguity (e.g. they change when Li is 
changed), while in their sum Mb = M^'^' +M^^' the ambiguity is reduced to 1 //n^. 

• In the present formulation of the effective theory, the 1/m-terms approach the continuum 
with an asymptotic rate oc a, in contrast to the leading order terms where this is [Q]. 

• Let us comment just on one numerical result at that point. The computation of a^"(Mi), 
eq. involves the difference of F^'"(2L) - F'^'"(L), where power divergent contributions 
cancel. As a typical case we choose L/a = 12 , T /a = 6, and the static action HYP2 (see [^), 
where our simulations yield a^rf "(2L) = 0.563 1 (6) , a^T\'''{L) = 0.5595(2), demonstrating 
a considerable cancellation. A detailed account of numerical results is presented in [^. 
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